Detecting phonon blockade with photons 
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Measuring the quantum dynamics of a mechanical system, when few phonons are involved, remains 
a challenge. We show that a superconducting microwave resonator linearly coupled to the mechanical 
mode constitutes a very powerful probe for this scope. This new coupling can be much stronger 
than the usual radiation pressure interaction by adjusting a gate voltage. We focus on the detection 
of phonon blockade, showing that it can be observed by measuring the statistics of the light in the 
cavity. The underlying reason is the formation of an entangled state between the two resonators. 
Our scheme realizes a phonotonic Josephson junction, giving rise to coherent oscillations between 
phonons and photons as well as a self-trapping regime for a coupling smaller than a critical value. The 
transition from the self-trapping to the oscillating regime is also induced dynamically by dissipation. 

PACS numbers: 85.85.+J, 03.67.Bg, 42.50.Ar, 03.75.Lm 



I. INTRODUCTION 

Since the early days of quantum mechanics the 
crossover to the classical world and the possibility that 
macroscopic objects could exhibit quantum behavior has 
attracted continuous interest. The tremendous progress 
in nano-fabrication capabilities have made these ques- 
tions amenable to experimental testing. Stimulated 
by the ideas of Leggett^, macroscopic quantum effects 
were first explored in Josephson junctions^ and nano- 
magnets^. In the recent past, the field of Nano Elec- 
troMechanical Systems (NEMS) has received much at- 
tention as a very promising ground for the investigation 
of these questions^^ and the observation of a mechani- 
cal oscillator in its quantum ground state^^^ is one of the 
most important achievements reached so far. 

The coupling of a quantum nano-mechanical oscillator 
to a qubit makes NEMS also suitable systems to explore 
the physics of circuit/cavity-QEEU^. Different schemes 
have been proposed including coupling to Cooper-pair 
boxes^^^^ and phase qubits^^^^^. Recently, coupling to 
a Cooper-pair box has been realized experimentally^^. 
Among the numerous interesting aspects of circuit-QED 
realized with mechanical resonators, here we want to ad- 
dress the phenomenon of phonon blockade^ which was 
considered recently in Ref. 20, extending to NEMS the 
original ideas put forward with photons in cavity- QED 
systems^^. The blockade effect arises because the co- 
herent coupling of the harmonic (photonic/phononic) 
mode with the (solid-state) atom leads to an effective 
non-linearity. For sufficiently strong coupling, the non- 
linearity is such that, upon external driving, the number 
of excitations of the oscillator never exceeds one. Obser- 
vation of phonon blockade in a nano-mechanical oscillator 
would be a clear evidence of its quantum nature. 

There is, however, an important challenge which needs 
to be tackled towards the demonstration of phonon block- 
ade: its detection. The motion of a mechanical oscilla- 
tor close to its ground state is tiny and strong ampli- 
fication of the transduced signals must be applied. On 



the other hand, amplification will inevitably add classi- 
cal correlations to the signal, thus disguising quantum 
correlations which are usually needed for demonstrating 
phonon blockade. 

In this article we show that these problems may be 
overcome if a Superconducting Microwave Resonator 
(SMR) is coupled linearly to the mechanical oscillator 
for transducing its motion into an electric signal. SMRs 
have proven to be nearly ideal quantum oscillators with 
easily tunable quality factors and they can be efficiently 
coupled to coherent quantum circuits^^^^. Our analysis 
is amenable to experimental verification. Recently, good 
coupling of a nano-mechanical resonator and a SMR has 
been demonstrated^^^^^^. Moreover, despite the difficul- 
ties in dealing with a microwave field at the level of a 
single photon, recently it has been shown that time cor- 
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FIG. 1. (a) Scheme of the system. A mechanical resonator 
is capacitively coupled to an artificial atom to induce a non- 
linearity and to a superconducting microwave resonator to 
detect phonon blockade, (b) Schematics of the energy spec- 
trum of the system for 77 ^ ^. The first excited states are the 
maximally entangled Bell states |^ ) and are off resonance 
with the states with more than one excitation. When the sys- 
tem is excited at the frequency ujd — uJo^g^ phonon blockade 
is observed through the photon statistics. 



relation functions for the cavity field can be measured ac- 
curately^— . Very recently, microwave photon blockade 
has been realized and measured experimentally^^. Our 
detection scheme takes fully advantage of the experimen- 
tal capabilities of the SMR to probe the quantum state 
of the mechanical resonator with a linear coupling. 

The hybrid phonon-photon system that we analyze in 
this work goes beyond a mere detection scheme for the 
phonon blockade. In the second part we present the prop- 
erties of this phonotonic junction. For a mechanical os- 
cillator initially driven in an exited state we show that 
two regimes exist, depending on the ratio between the 
anharmonicity and the coupling compared to a critical 
value. For a small coupling, the phonons are trapped 
in the mechanical oscillator whereas coherent oscillations 
between the phonons and the photons appear when in- 
creasing the coupling, leading to quantum revivals in the 
photon statistics. We show that this transition can be in- 
duced by the dissipation and take place during the tem- 
poral evolution of the system. 

The paper is organized as follows. In the next Sec- 
tion we define the model for the system depicted in Fig(T] 
and introduce the equations governing its dynamics in 
the presence of damping. In Section IIIII we discuss how 
to detect phonon blockade by means of a measurement 
of the photon statistics. In Section |IVl we present the 
properties of this phonotonic junction with the transi- 
tion between a self-trapping regime and a regime with 
coherent oscillations. We summarize our results in the 
concluding Section. 



II. THE MODEL 

The system we have in mind is depicted in Fig. [TJ It 
is a mechanical resonator coupled capacitively to an ar- 
tificial atom and a SMR. The non-linearity, leading to 
phonon blockade, is induced by coupling the oscillator to 
a superconducting qubit. In the case of a Cooper-pair 
box, one gets the Jaynes-Cummings model. A Cooper- 
pair box molecule^^ may be considered as well in order to 
increase the non-linearity. In the present work, however, 
we do not deal with a specific choice of the superconduct- 
ing nano-circuit; the important ingredient is the gener- 
ation of the Kerr Hamiltonian proportional to a^a^aa^"^, 
where a and a^ are the phonon annihilation and creation 
operators of the mechanical oscillators. After the adia- 
batic elimination of the qubit, the Non-linear Mechanical 
Resonator (NMR) is described by the following effective 
Hamiltonian 



^NMR = hwrO/ a + hrfo/ a 



aa. 



(1) 



The mechanical resonator is supposed to be in the quan- 
tum regime (the bare frequency is in the GHz range). 
For a strong coupling of hundreds of MHz, the Kerr non- 
linearity strength 77 is of the order of MHz. The coupling 
between the NMR and the SMR is obtained from circuit 



theory. The cavity is modeled by an array of LC cir- 
cuits^^. The Hamiltonian of the fundamental mode of 
the SMR is i^sMR = fujjjj^h^ where ujc/^n is the mode 
frequency and h {l)^) the corresponding photon annihila- 
tion (creation) operator. 

The mechanical resonator is kept at a fixed potential 
Vg with respect to ground and the coupling is realized 
through a localized capacitance Cg = Cg -\- {a -\- a^)Cg, 
resulting in two coupling terms in the Hamiltonian: a 
radiation pressure term CgV^^^ {a -i- a^)b^b and a linear 

term CgVgVr,m.s.{^~^^^){^~^^^)^ coming from the electro- 
static energy of the coupling capacitor (K.m.s. is the root 
mean square of the zero-point voltage fluctuation of the 
SMR at the position of the coupling capacitance). These 
two terms can have very different orders of magnitude 
for typical values of Vg and K.m.s.^^- We focus on the 
case where the gate voltage is much larger than the root 
mean square of the voltage quantum fluctuations inside 
the SMR, and thus neglect radiation pressure. In the 
rotating wave approximation, the coupling reads 



Hint = % (a^b^ab^) , 



(2) 



where g = CgVgVr.m.s./^- For realistic parameters, the 
coupling may be of the order of MHz. The mechanical 
resonator is driven by a weak oscillating current, at the 
frequency ujd/27r^ in the presence of a static magnetic field 
perpendicular to the plane of the circuit. This driving is 
modeled by the Hamiltonian 



H(^^^ = he [a^e *' 



iujdt 



ae 



iuJdt\ 



(3) 



where the strength e is proportional to the current am- 
plitude and the magnetic field. Driving the mechanical 
oscillator can in principle act as a direct driving also for 
the SMR through the capacitive and inductive couplings. 
However, the resulting driving is negligible. Indeed, on 
the one hand because the displacement is tiny the electro- 
motive force induced in the SMR is well below K.m.s. and 
on the other hand for spatially separated resonators the 
mutual inductance is suppressed. The total Hamiltonian 
of the system reads H = i^NMR + ^smr + ^int + ^drv 
We choose the working point at Ur = ujc = 000- The 
mode frequency of the cavity can be tuned by adding a 
SQUID at the end of one arm to change the boundary 
condition^^. In the rotating frame of the driving, the 
total Hamiltonian reads 

H = h{uJo — ujd){ci ci -\- ¥b) -\- hrfo/ 0/ aa 

+ hg{a^b + ab^) + he{a^ + a). (4) 

A linear coupling between driven non-linear oscillators 
can also be obtained in cavity-QED^^. 

The finite lifetime of the phonons and the photons is 
taken into account through the Lindblad operators L^ 
and Lc of the resonator and the cavity, respectively^^ (p 
is the density matrix of the whole system) 

Lrp = ^7r {2apa^ — a^ap — pa^a) , (5a) 

Lcp = he (2^/>^^ - b^bp - pb^b) , (5b) 



where the damping rates ^r,c = ^r,c/Qr,c are the inverse 
of the phonon and photon Hfetime and are defined by the 
quahty factors Qr,c- The dynamics of the system is then 
governed by the master equation 



dtp{t) 



1 
ih 



[H,p{t)]+Lp{t), 



(6) 



where L = L^ + Lc is the total Lindbladian. 

In principle, some noise is introduced by the voltage 
source used to keep the mechanical resonator at Vg . How- 
ever, following Ref. 39, the dominant Lindblad operator 
corresponding to this source of noise for the case of a 
Markovian environment is found to be of the order of 
r]gRCgVr.m.s./Vg {R is the internal resistance of the volt- 
age source), which is negligible due to the very small 
value of gRCg. In the case of slow voltage fluctuations, 
one can also neglect noise effects, as we discuss below. 



III. DETECTION OF PHONON BLOCKADE 

Blockade is possible only if the non-linearity of the en- 
ergy spectrum is larger than the state linewidth, namely 
V^9 > lr,c' This condition imposes the quality factors 
to be at least several thousand, which is within the ex- 
perimental capabilities. Throughout the present work we 
chose the following parameters for the numerical simula- 
tions: u;o/27r = 1 GHz, r]/27r = 10 MHz, g/27r = 1 MHz, 
Qr = 10^, and Qc = 10^. Fig. [3] shows that the results 
are still valid for quality factors (5r,c = 3000, when the 
damping becomes of the order of the coupling. Moreover, 
the NMR is driven at the resonance ujd = ojq — g with the 
amplitude e/27r = 0.1 MHz. By analyzing the time traces 
of the photon and phonon populations (not shown here) 
it is possible to note that the cavity closely resembles the 
dynamics of the NMR, with equal steady state phonon 
and photon numbers close to 0.25. The population of 
the state with two phonons or two photons is strongly 
suppressed, implying both phonon and photon blockade. 
In the rest of the paper, we will show by solving Eq. (|6|), 
that by means of the detection of photon correlations it 
is possible to extract unique information on the phonon 
statistics. 

In order to ascertain the accuracy of the proposed de- 
tection scheme, it is not sufficient to see a correlation in 
the populations. We analyze the statistics of the excita- 
tions by means of the second order correlation function^^ 



-,(2) 



(r) 



lim 



{yKt)y{t))' 



(7) 



.(2)/ 



where y = a or b. The value of g^lir) is comprised 
between and 2 and tends towards unity for long time 
difference, where the coherence is lost. The value at co- 
inciding times r = reffects the statistics of the field: 
a value of ga (0) <C 1 corresponds to antibunching, and 
serves as the signature of phonon blockade. 
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FIG. 2. (Color online) Spectroscopy of the energy spectrum. 
Steady state phonon number (solid line) and photon number 
(dashed line) as a function of the driving frequency. The 
different peaks correspond to the excitation of the states with 
one phonon |l,nc). Inset: Second order correlation functions 
at coinciding times for the phonons Qa (0) and for the photons 
g^ (0) . Anti-bunching occurs at cj^ = cjq ± S' where phonon 
and photon blockade take place. Photon bunching occurs 
when the states with many photons are excited. 



To understand how to induce and detect phonon block- 
ade, we look at the energy spectrum of the undriven 
Hamiltonian H in the Fock basis |nr,nc), where rir^c 
is the phonon and photon number respectively. The 
total number of excitation a^a + b^b being conserved, 
the spectrum can be decomposed on the subspaces de- 
fined by a given number n = Ur -\- Uc of excitations 
{|/c,n — /c), k = 0, ...,n}. For one excitation n = 1, 
the eigenstates are the maximally entangled Bell states 



|vI/±) = -i^(|0,l>±|l,0», 



(8) 



with the energy hcuo ± hg. If there is a small detuning 
between the two resonators, the eigenstates are rotated 
by an angle {ur — ujc)/^g- For higher excitation numbers 
the ladder structure depends on the ratio rj/ g between the 
anharmonicity and the coupling. In the limit of a strong 
non-linearity r]^ g^ the spectrum is composed of two en- 
tangled states |0, n)±|l, n — 1) at uj{)^^/ng and n — 1 pure 
states |2, n — 2), . . . , |n, 0) located at flcjo + ^^(^^ — 1)^, 
see Fig. [1] (b). This non-linear spectrum allows for ex- 
citation blockade, since the energy of the state |^^) is 
not resonant with higher states. If the system is excited 
at the frequency uod = ^q ^ g-, only one excitation is cre- 
ated, symmetrically shared between the NMR and the 
SMR. The eigenstate being the maximally entangled Bell 
state 1^^), the photons have the same dynamics as the 
phonons and the cavity constitutes consequently a reli- 
able measurement device to detect the state of the res- 
onator through the photon statistics (see Fig. H]). 

The energy spectrum can be probed with the response 
of the system to the driving current when the driving 
frequency is tuned, as depicted in Fig. [2l In order to 
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FIG. 3. (Color online) Second order correlation functions 
against the quality factors. The phonotonic Josephson junc- 



tion works for quality factors down to Qr 



W, where the 



damping rates are comparable to the coupling strength. 

excite states with one phonon such as |1, n — 1), the driv- 
ing frequency is fixed to uod = o^o ^ ^/v^- These val- 
ues correspond to the peaks in the excitation numbers 
of Fig. [21 Compared with the dependence of the g)^ ^(0) 
on ujd (see the inset), it shows that blockade occurs at 
<^d — ^0 = '^9 where the second order correlation func- 
tion is minimized. This minimum protects the blockade 
phenomenon against slow fluctuations of the gate voltage 
Vg^ or equivalently the coupling g. Indeed, if the driving 
is kept at a fixed frequency U0d = oo^^g and one looks at 

(2) 

the dependence of g)^ ^(0) on a coupling g(t) fluctuating 
around the value ^, one obtains a behavior analogous to 
the inset of Fig. [2l namely a minimum at g = g. The 
averaged value of the second order correlation function 
over a Gaussian distribution of coupling strengths is not 
affected and both phonon and photon blockade is thus 
insensitive to slow gate voltage fluctuations at first or- 
der. For higher values of n, the photon number increases 
and the second order correlation function of the cavity 

(2) 

gl ^(0) tends to 2, indicating photon bunching. 

When the anharmonicity is very large, the NMR be- 
haves like a two-level system and can be described by re- 
placing the operators a and a^ by the the ladder operators 
a- = |0)(1| and cr+ = |1)(0|, respectively. The Hamil- 
tonian is then reduced to an effective Jaynes-Cummings 
model, Hjc = hg{a-b^ + <^+^)- ^^ lowest order in g/r]^ 
the asymptotic expression of the second order correlation 
functions is given by 



g'a 



(^)(0) = (46/7^)^ 



9b 



i'\0) = {2e/g)\l^4g/rj). (9) 



The comparison with the numerical results is presented 
in Fig. m where the correlation functions are plotted as a 
function of the anharmonicity for different values of the 
coupling. For a sufficiently large anharmonicity the prop- 
erties of the SMR become 77-independent with a strong re- 
duction of the phonon correlation function ga (0). Pho- H = huJo{a^a -\- b^b) + hrja^a^aa + hg{a^b + ab^). This 



FIG. 4. (Color online) Steady state phonon and photon num- 
bers (a) and second order correlation functions (b) as a func- 
tion of the Kerr non-linearity. The Bell state |^~) ensures 
a perfect match between the phonon and photon numbers. 
Many excitations are generated when the non-linearity is com- 
parable with the damping rates of the resonators, where the 
excitations are not anti-bunched anymore. The correlators 
are plotted for different values of the coupling, from top to 
bottom, g/27T = 1, 2, 5, 20 MHz respectively. The correla- 
tion of the NMR is essentially coupling-independent while for 
large 77 the correlation of the SMR saturates to a constant 
value (2e/g)^. The corresponding solutions Eq. (|9|) in the 
limit 7] ^ g are plotted in dotted lines. 



ton blockade is enhanced when the coupling increases. In 
the opposite limit of a small anharmonicity, the driving 
generates many excitations. The transition from photon 
anti-bunching to photon bunching is also observed when, 
for a fixed anharmonicity, the coupling decreases or the 
driving amplitude increases. The former is due to the 
degeneracy of the Bell states for small coupling g ^ jr,c 
and the latter is because of the effective level broadening 
due to the driving. 



IV. THE PHONOTONIC JUNCTION 

For an isolated mechanical oscillator {rj = g = jr = ^) 
initially in its ground state |0)r, the resonant driving 
{ud = cjOr) generates a coherent state \a = —iet)r = 
ex.p[—iet{a^ -\-a)]\0)r after a time t (see Fig. [5] (a)). In the 
presence of dissipation, at zero temperature, the steady 
state is the coherent state \a = —i2e/jr)r-^^' The synthe- 
sized coherent state |a)r, with a tunable phonon num- 
ber |ap, will be used as the initial state in the follow- 
ing. Once the phononic state is prepared, the driving is 
turned off while the coupling to the qubit and the SMR 
are switched on. The latter can be performed by tuning 
the qubit frequency closer to the resonance and putting 
the gate voltage Vg on, respectively. The system is then 
governed by a two-sites Bose-Hubbard like Hamiltonian 
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FIG. 5. (Color online) (a) Population distribution of the 
mechanical oscillator after being driven to a state with 10 
phonons (points). The statistics is compared to a coherent 
state with the same mean value (line), (b) Average value 
of the imbalance z as a function of the parameter A from 
the classical dynamics without dissipation, starting from the 
initial state z{0) — 1. The symmetric junction corresponds 
to interacting photons with a strength 77. (c) Time evolu- 
tion of the population imbalance in the self-trapping regime 
(?7/27r = 10 MHz and g/2T^ — IMHz), the oscillating regime 
\r]/2n = IMHz and g/2n = 10 MHz). The dynamical tran- 
sition is observed with a symmetric junction, where the pho- 
tons also interact with a strength ?7/27r = 10 MHz. The cou- 
pling is equal to g/2Ti = 10 MHz and the quality factors are 

Or.c = 10^ 



is similar to the Bosonic Josephson junction, realized 
with a cloud of cold atoms in a double well potential^^ . 
Since only one of the two resonators is non-linear, we 
are dealing with an asymmetric junction. The symmet- 
ric case, where the photons are also interacting, can be 
obtained by adding a qubit in the SMR to generate a 
non-linearity^^. To present the properties of this junc- 
tion, we start with a classical description in the absence 
of dissipation. 

The dynamics of the phonotonic junction can be de- 
scribed in terms of the imbalance z = {ur — Uc) / {nr^ric) 
between the phonons and the photons and the phase dif- 
ference Lp = diTg{a^b) between the two resonators. In 
the absence of dissipation, the total number of particles 
{rir + Tic) is conserved, equal to |ap. The classical equa- 
tions of motion for (n^), (^c), and (a^b) can be expressed 
in terms of z and (p. The system is then classically gov- 
erned by the following set of nonlinear differential equa- 
tions 



z{r) = y/l — z'^{r) sin(/9(r), 

z{t) 



<p{t) = A[l + z{t)] 



n/1 - zHr) 



cos(/?(r). 



(10) 

(11) 



large couplings, coherent oscillations take place between 
the phonons and photons. When the interaction between 
phonons is larger than the coupling, the oscillations are 
frozen and the particles are self-trapped. The transition 
between these two regimes is controlled by the critical 
parameter Ac = 2, corresponding to a critical coupling 



c{t) = \r]{nr^nc){t). 



(12) 



where the time has been rescaled to r = 2gt. The pa- 
rameter A = ri{nr + Tic) /2g leads to two regimes^^. For 



These two regimes are presented in Fig. [5] (b). The 
time average of the imbalance vanishes in the oscillating 
regime A < Ac and tends to unity in the self-trapping 
regime A > Ac. The imbalance is also presented for a 
symmetric junction, where the photons are also interact- 
ing with a strength rj. For a symmetric junction, the 
imbalance follows the dynamics of Eq. (p!Q|) while con- 
cerning the phase the first term in the right-hand side of 
Eq. (pTj) has to be replaced by 2Kz{t). This gives rise to 
a sharp transition at A = 1. 

In the presence of dissipation, the critical coupling de- 
creases with the characteristic rate 7^^c- If the coupling 
is initially larger than the critical coupling gc{^)^ the 
system starts in the self-trapping regime. After a time 
to ~ ln[^c(0)/^]/7r,c7 the coupling becomes larger than 
the critical coupling gdto) and coherent oscillations take 
place between phonons and photons. Quantum revivals 
can then be seen in the cavity. This dynamical tran- 
sition induced by dissipation is sharper when the pho- 
tons are also interacting^^. The transition is presented 
in Fig. [5] (c) for a symmetric junction. 



V. CONCLUSION 

In conclusion, we have shown that coupling a SMR 
to a NMR is a powerful tool to detect phonon blockade 
and generate entanglement between phonons and pho- 
tons. The main reason for the accurate detection when 
few phonons are involved is the formation of Bell states 
between the two resonators, ensuring a perfect match 
between the phonon dynamics and the photon statistics. 
The phonotonic Josephson junction takes advantage of 
the recent experimental capabilities with microwave pho- 
tons. The simulations, obtained in the framework of the 
quantum master equations, demonstrate that our pro- 
posal is compatible with the current experimental capa- 
bilities. Our new detection scheme realizes a phonotonic 
Josephson junction and its applications go beyond the 
blockade regime. A rich physics stems from this device 
which can be used for instance to observe phonon las- 
ing through microwave photon lasing or the dynamical 
Casimir effect if the gate voltage is strongly modulated. 
The phonotonic Josephson junction constitutes a build- 
ing block towards the use of NEMS as quantum buse&^. 
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